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Abstract. We propose a multi-scale approach to understanding physics related to the ion/proton-beam cancer therapy and 
vQ ' calculation of the probability of the DNA damage as a result of irradiation of patients with energetic (up to 430 MeV/u) ions. 

This approach is inclusive with respect to different scales starting from the long scale defined by the ion stopping followed by a 
smaller scale defined by secondary electrons and radicals ending with the shortest scale defined by interactions of secondaries 
with the DNA. We present calculations of the probabilities of single and double strand breaks of the DNA and suggest a way 
r^ ' of further elaboration of such calculations. 
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t>-»' Ion-beam cancer therapies are being used more and more as favorable alternatives to the conventional photon therapy, 

1-5 i also known as radiotherapy [ 1 ] ' . Their advantages related to at the fundamental difference in the linear energy transfer 

(LET) by a massive projectile as compared with massless photons, namely by the Bragg peak depth-dose distribution 
for the ion. It is due to this peak, can the effect of irradiation on the tissue be more localized increasing the efficiency 
of treatment and reducing the side effects. In order to plan a treatment, a number of physical parameters, such as the 
^ \ energy of projectiles, intensity of the beam, time of exposure, etc., ought to be defined. At present their definition 

is based on a set of empirical data and experience of personal. Moreover, the optimization of treatment planning 
requires understanding of microscopic phenomena, which take place on time scales ranged from 10^^^ s to minutes or 

^k^ . even longer times. Many of these processes are not sufficiently studied. Thus, a reconstruction of the whole scenario 

explaining, qualitatively and quantitatively, the leading effects on each structural level scale presents a formidable task 
not only for physics but also for chemistry, biology and medicine. 

^Q ' The ultimate effect of the beam therapy is due to the DNA destruction and subsequent killing of cells as a result 

f^*) \ of energy deposition by the projectiles [1]. Most of energy deposition by the ion is due to ionization of a medium it 

traverse through, which by about 75% consists of liquid water The secondary electrons formed in this process are 
believed to be mostly responsible for the DNA damage, either by directly breaking the DNA strands, or by reacting 
with water molecules producing more secondary electrons and of free radicals, which can also damage the DNA. One 
^ \ can also speculate about heating of the medium in the ion tracks making the DNA more vulnerable to damage, if 

not melting it. Among the DNA damage types, we emphasize single strand breaks (SSB's) and double strand breaks 
(DSB's). The latter ones are especially important because they represent irreparable damage to the DNA. After a fast 
ion beam enters the tissue, many processes take place on different temporal and spatial scales until tumor cells die. 
The goal of our approach is to analyze these processes and select the main physical effects which are responsible for 
the success of the ion-beam therapy. It turns out that several aspects play an important role as is illustrated in Fig. 1 . 

The driving force of the ion-beam therapy is a propagation and stopping of the incident ions in the tissue. Depending 
on the initial kinetic energy, the ions penetrate to a certain depth and produce a Bragg peak at the end of their 
range. Many works were devoted to this problem including both deterministic and Monte Carlo methods, see e.g. 
Ref.[2]. Using the information about cross sections of atomic (such as ionization of water molecules) and nuclear 
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' We limit the references to either reviews or the most recent works that have caught our attention and referred to other literature on the subject. 
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FIGURE 1. Schematics of the multi-scale approach. 

(such as nuclear fragmentation of projectiles) processes as an input, these models give very good predictions of 
all characteristics of the Bragg peak, its position, height, tail, etc. These models provide a reasonable information 
on the energy deposition on a mesoscopic scale of about 0.1mm, which is sufficient for the treatment planning. 
On the other hand, this information is not sufficient to analyze the processes taking place on a microscopic scale. 
The energy of ions changes from the initial energy in the range of 200-430 MeV/u down to about 50 keV/u. The 
next scale is defined by the secondary electrons, which are produced as a result of ionization of molecules of the 
medium and by radicals also formed as a result of energy loss by the projectile. The maximum energy on this scale 
hardly exceeds 100 eV and the displacement is of the order of 10-15 nm. So far, we have only considered ionization 
of water, a surrogate for a biological medium, as a single process taking place on this scale, leaving other atomic 
interactions for later consideration. Ionization of the medium is the leading process in accounting for the energy 
deposition by the ions we believe that the secondary electrons are mainly responsible for the DNA damage. Interaction 
of electrons and radicals with DNA also happens on a nanometer scale, and many works are devoted to study these 
interactions [10, 11, 12, 13, 14, 15]. 

The main event of this scale is a diffusion of free electrons and radicals in the medium. Many chemical reactions 
take place as well. They are also important for estimates of the DNA damage since they define the agents interacting 
with the DNA. Again, this aspect attracted plenty of attention of Monte Carlo simulations adepts [5] who using various 
SDCS for ion and electron energy loss (including the effects of the medium [6, 7, 8, 9]) trace the electrons and other 
species through the medium up to their interaction with the DNA. In this paper, we present an approach to calculations 
that can be done on this scale without using Monte Carlo simulations. 

. The main input for our approach is the singly differentiated cross section (SDCS) of ionization of water [3, 4]. In 
this work, we use the experimental results of Ref. [10, 11, 12]. 

PHENOMENA-DEFINED SCALES 
Ion stopping and production of secondary electrons 

The energy spectrum of secondary electrons was analyzed in our previous works, [3, 4]. It plays a major role in the 
energy loss by projectiles and therefore determines to a large extent all characteristics of the Bragg peak. As follows 
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from Eq. (5), it is also important for the DNA damage calculations (see below). 

Physically, the SDCS is determined by the properties of the medium, and since we use liquid water as a substitute for 
biological tissue, it is determined by the properties of water molecule and the properties of liquid water as continuous 
medium. All this information is contained in real and imaginary parts of electric susceptibility of liquid water. This 
approach can be generalized for any real tissue. In order to do that, the same quantities, such as SDCS, of this medium 
have to be used. 

In Refs. [3, 4], we used a semi-empirical parametrization by Rudd [16] for the SDCS and obtained the position of 
the Bragg peak with a less than 3% discrepancy as compared to Monte Carlo simulations and experimental data [3]. 
However, contrary to the calculations of total cross sections, or the Linear Energy Transfer (LET), where there is 
integration in W, the calculations of the DNA damage may be more sensitive to the shape of the SDCS at small 
energies, which for the liquid water is significantly different from that of water vapor [17, 18, 19, 20]. 

The SDCS is a function of velocity of the projectile and, since the ions are quite fast in the beginning of their 
trajectory, it has to be treated relativistically. In Ref. [3], this issue has been solved by "relativization" of the Rudd 
parametrization by fitting it to correct Bethe asymptotic behavior at the relativistic limit. 

Another important issue related to SDCS is the effect of charge transfer that is due to picking-up electrons by the 
initially fully stripped ions (such as '^C^+) as they slow down in the medium. Since the SDCS is proportional to the 
square of ion charge, its reduction strongly influences such characteristics as the height of the Bragg peak, secondary 
electron abundance, etc. In Ref. [3], we solved this problem by introducing an effective charge taken from [21]. As a 
result, the effective charge of the ^^C^+ near the Bragg peak is about +3 rather than +6. 

Even after the relativistic treatment of the projectile and the introduction of effective charge, the profile of the 
Bragg peak obtained in our calculations was substantially higher and narrower than those obtained by Monte Carlo 
simulations or experimentally. The main reason for the discrepancy was that our calculations were performed for a 
single unscattered ion, while in simulations as well as in experiments the ultimate results are a combination of many 
ion tracks with a significant spread in energy and position due to multiple scattering by water molecules. After we 
took into account straggling of the ions, the shape of our Bragg peak matched the shape predicted by the Monte 
Carlo simulations with nuclear fragmentation channels blocked (see for details in this collection the contribution of 
Ref. [22]). 

The nuclear fragmentation in the case of carbon ions is quite substantial and should not be neglected. In principle, we 
can include the beam attenuation due to nuclear reactions given the energy dependent cross sections of these reactions. 
Then we would be able to reproduce the attenuation of the ion beam, secondary electron production due to different 
species, the spread of the Bragg peak due to different penetration depths of different species, and the tail following the 
Bragg peak due to light products such as protons and neutrons. All these complications, however, were beyond our 
primary goal of gathering most significant effects together and we leave it to the future refining projects. We should 
mention that a successful treatment of nuclear processes has been done by the Monte Carlo simulations [2]. 

Another process not accounted for in Ref. [3] is an excitation of water molecules by the ions. This effect contributes 
in the energy loss by the projectiles and shifts the Bragg peak towards the source. Even though no secondary electrons 
are produced in this case, the excitation channel may be important for the DNA damage since excited water molecules 
may dissociate producing free radicals that, in their turn, may damage the DNA. 

Thus, the ionization energy loss by the ions in liquid water is the dominating process for the ion stopping and the 
energy spectrum of the secondary electrons. Additional energy losses are associated with excitation of water molecules 
leading to the production of free radicals. The SDCS defines both the shortest scale related to the ion propagation and 
provides the initial conditions for the next scale related to the secondaries. 



Propagation of the secondary electrons 

Even though the SDCS that we have used in [3] may not be quite adequate, this distribution gives some important 
predictions that agree quite well with other calculations and measurements. Indeed, the average energy of the secondary 
electrons, 

I r da(W,T) 
{W) = — W \' ' dW, (1) 

c7t Jo dW 

in the vicinity of the Bragg peak (T « 0.3 MeV/u) is about 45 eV. This value limits possible further scenarios for such 
electrons. For instance, it was shown that such electrons may excite or ionize another water molecule, but, most likely, 
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only once, and the next generation of secondary electrons is hardly capable of ionizing water molecules [3, 4]. This 
puts a cap on the amount of produced secondary electrons. 

The secondary electrons propagate in the same medium as the ion, and interaction with the medium is again 
determined by the SDCS with electrons being projectiles. The interaction can be elastic or inelastic, and there is a 
probability of stumbling on a DNA and causing its damage. 

The angular distribution of the secondary electrons at energies about and below 45 eV are rather flat [23]. Therefore, 
to the first approximation, we can consider a Brownian motion of secondary electrons and use random walk or simply 
diffusion to describe their propagation through the medium from a point on the ion trajectory where they became 
unbound. The probability density to diffuse through the distance r after k steps is given by the following [24], 

In this equation, the mean free path I is the average distance that the electron passes between two consecutive elastic 
collisions. The elastic mean free path is determined by the SDCS of electrons and we use the results of Ref. [8]. 
Its typical values are 0.3-1 nm for electron energies around 30 eV. The mean free path for inelastic collisions /,„ 
is typically many (about 10-20) times longer So, we assume that electrons mainly experience elastic collisions and 
inelastic processes are included via attenuation factor exp(— ZA:/Z,„), accounting a distance wandered kl . 

Both elastic and inelastic mean free paths depend on the energy of the wandering electron. This energy is changing 
gradually, and strictly speaking, the mean free path is a function of the initial energy and the number of steps. However, 
the estimates of ref.[7] suggest that the electron energy loss over typical distances of 6-10 nm is not very large. 
Therefore, we will use some average energies and assume that they do not change during the diffusion. 

Evaluation of the DNA damage 

A DNA damage, such as a SSB, is a result of a sequence of mutually independent events. First, secondary electrons 
with a certain kinetic energy W are produced at a certain depth x. Later on they interact with a surrounding medium 
via elastic and inelastic collisions, and gradually lose energy until become bound. Depending on the electron's energy, 
momentum and position, there is a chance that an electron stumbles on a DNA molecule and damages it. Fig. 2. 
To evaluate the DNA damage by secondary electrons we propose the following model. First, we represent the DNA 
molecule by a single convolution, i.e. a cylinder of length 3.4 nm and radius 1.1 nm (these parameters are well 
estabilished experimentally). The reason for that is that the DSB's are defined as simultaneous breaks of both DNA 
strands located within the single convolution. Second we assume that this cylinder is irradiated by the flux of secondary 
electrons produced by the ion traversing the medium at a certain distance from the DNA. This geometrical picture is 
schematically shown in Fig. 2. It is clear that the mutual location and orientation of the DNA cylinders and ion 
trajectories are randomly distributed. Therefore, the total damage to the DNA molecules can be calculated by averaging 
over all the possible configurations. From that scheme stems that any distance r between a point on the track and a 
point on the cylinder is given by 

r^ = (flcos^-p)^ + (flsin^- ^sina)^ + (z-zq- Ccosa)^ (3) 

The flux of secondary electrons through a unit area at a distance r from the production point is 

d,(r,W,W/) -DV/?(fc,r)^^(C,W)AW , (4) 

where D = kfi /6 is the diffusion coefficient multiplied by the average time of wandering. The last factor in Eq. 4 is 
the number of electrons with energies between W and W + AW produced from the segment of ion trajectory dC,. It 
can be obtained from the singly differentiated ionization cross section (SDCS), as explained in refs. [3, 4]. Wf is the 
electron's energy when it reaches the surface after traveling from the track. Finally, we assume that the number of the 
SSB's within the DNA cylinder is proportional to the number of electrons crossing the surface, 

NssB = TssB {Wf)Y^j d^dA ■ *^ (r, W,Wf) , (5) 

k •' 
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FIGURE 2. Geometry of the model: z is the cyUndrical axis of the DNA convolution and x is chosen to be parallel to PQ, the line 
of closest approach between z and the beam (ortogonal to both), of length p, at distance Z{) from the center of the convolution O. ^ 
is the coordinate of any point in the beam with respect to P and a is the angle between the beam and z- In the right panel we show 
a projection on the xy plane, where a = l.lnm is the DNA radius and <j> is the polar angle defining a point on its surface. 

where the integration is done over the surface of cylinder and the ion trajectory. The unknown quantity TsssiWf), 
that for the moment we assume as a constant, should be determined from experimental data. The calculation for a 
general case as shown in Fig. 2 is rather cumbersome. Therefore, let us set zo to zero, and consider separately two 
limiting cases, the "parallel" case when a — 0, and the "normal" case when a = k/2. In the parallel case, the cyUnder 
containing the DNA convolution is parallel to the ion track and p is the distance between the axis of the cylinder and 
the track. In the normal case, the axis of this cylinder is perpendicular to the ion track and when zq ~Q, p is again the 
distance between the axis of the cylinder and the track and the beam projects along p to the center of the cylinder. In 
both of these cases we need to set the Umits for angular integration over (p. 

Looking from any point of the ion track, there are two surfaces of the DNA cylinder: the "front" or "face" surface 
and the "back" surface (see Fig. 2). In our model, if a wandering electron hits the face or the back surface, it may cause 
a strand break with a certain probability. Therefore we simply add the probability of a SSB due to electrons striking the 
back surface to that for the face surface regardless of directions of their motion leaving an introduction of attenuation 
mechanism accounting for electron passage "through the DNA" to an extension of this model. 

The results of the integration are shown in Fig. 3. All curves show the dependence on the distance p from the DNA 
to the beam. When the distance is large enough (> 3 nm) the normal and parallel cases coincide for the face side as 
well as for the back side, only at small distances there is a difference. This difference is significant only for back and 
face sides taken separately, but not for their combination shown in Fig. 4a. This means that the geometrical details 
of orientation of DNA segments with respect to the beam may not be so significant, since all varieties lie somewhere 
in the shadowed region in Fig. 4a between the two curves. For an analysis of a more general picture, some average 
curve (lying between the two curves in Fig. 4) should be used with p^ replaced by p^+Zq- In this calculations we took 
the value TssaiW/) ~ 5 x 10^'*, extracted from experimental data of Refs. [10, 11, 12], where SSB's and DSB's were 
induced by the electron beam of energy 0.1-30 eV. The density of the beam was such that only a small fraction of 
DNA molecules was irradiated. The numbers of SSB's caused by the secondary electrons depending on the distance p 
and the energy of the secondary electrons for the parallel case are shown in Fig. 5. The decline of the number of SSB's 
with increasing p is an explicit consequence of Eq. (2) with account of (3). The energy dependence is mainly caused 
by the dependence of mean free paths on the energy and the attenuation factor given above. This factor is heuristic and 
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FIGURE 3. A comparison of numbers of SSB's for parallel (dashed line) and normal (solid line) configurations on the face (a) 
and on the back side (b) for 20-eV electrons. 




FIGURE 4. A comparison of dependencies of overall (due to the whole surface of the cylinder) SSB's (a) and DSB's due to 
separate electrons (b) on distances to the DNA convolution in the parallel (dashed line) and normal (solid line) cases for 20-eV 
electrons. 



may have to be corrected later when the corresponding experimental or computational data are available. 

What about the double strand breaks? From Ref.[10, 11, 12], it follows that the DSB's produced by the electrons 
with energies over about 5 eV are caused in one hit, i.e., if a particular electron with a probability of about 0.0005 
causes a SSB, the same electron causes a DSB with a probability of about 0.0001. This is why the analysis of the 
probabilities of SSB's are so important. If the energy of the secondaries are high enough they give the probability of 
the DSB's after being divided by some factor of about 5. Therefore, Fig. 5 also gives the shape of the main contribution 
to the DSB's. This, of course, strongly depend on the electron density and energy. Many authors e.g. [8] use a criterion 
that electrons with energy above a certain threshold (say 30 eV) are needed to produce a DSB. At energies lower than 
5 eV the situation changes: one electron does not cause two breaks. Therefore, we need to calculate the number of 
DSB's caused by two different electrons. From the geometry of a DNA convolution we infer that the probability of a 
DSB is proportional to ^{N^sb face '^^ssb back) + j^ssB.faceNssBjmck- The numbers for the DSB's caused by different 
electrons in parallel and normal cases are shown in Fig. 4b. As well as in the case of SSB's the difference due to 
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FIGURE 5. The number of SSB's dependent on the distance to the DNA convolution, p, and the energy of secondary electrons 
(parallel case). 



geometry is not very significant. Even though the numbers of DSB's are many times smaller than those in Fig. 4a, 
this effect may be significant if the density of secondary electrons is large enough. According to our estimates in 
Refs. [3, 4], the density of the secondary electrons in the vicinity of the Bragg peak at therapeutic conditions is by 
about 16 orders of magnitude higher than the electron density in experiments of Ref. [10, 11, 12]. Therefore, this effect 
may be an important correction to the paradigm. 

This concludes our approach to calculations of DSB's and SSB's due to secondary electrons produced by the ions. 
What about secondaries produced by the electrons and free radicals produced by the ions? Can they be treated the 
same way as we have treated the secondary electrons in the previous subsections? 



Other secondaries 

Secondaries that can be treated almost in the same way as the secondary electrons that we just went through are 
OH- radicals. They are formed as a result of ionization of water molecules by the ion after dissociation of water ion 
into OH- and H+. These radicals are formed almost at the same place as the secondary electrons. The difference is, of 
course, in a different diffusion coefficient, and different time of getting to the DNA, which is by about 100 times longer 
than that for secondary electrons. Then, the DNA damage caused by OH- may also be different [25, 26]. Nonetheless, 
if the effect produced by OH- is an important player, this is a recipe of its inclusion. The same can be said about those 
free radicals that are formed as a result of excitation of water molecules by the ions. These radicals (OH- and H-) are 
also produced on the ion trajectory and can be treated in the same way as secondary electrons. 

The other secondaries, such as the second generation of electrons produced by the first generation via ionization of 
water, radicals produced as a result of this process and the radicals H- produced via interaction of secondary electrons 
with water molecules {e.g., through dissociative attachment) can be treated in the following way. Let the interaction 
that produces a "desired agent" happens at some point r'. Then the previous procedure has to be divided into three 
parts: diffusion of the secondary electron from the point of origin (the ion's trajectory) to F', an interaction that leads 
to the production of the agent at r', and the diffusion of the agent to the DNA cylinder. Then, perhaps cumbersome, 
integration over r' has to be performed. 

CONCLUSIONS 

Thus we presented a multi-scale inclusive approach to the physics relevant to the ion-beam cancer therapy. We intend 
to present a clear physical picture of the events starting from an ion entering a tissue leading to the DNA damage as a 
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result of this incidence. We view this scenario as a palette of different phenomena that happen on different time, energy, 
and distance scales. From this palette, we choose major effects that adequately describe the leading scenario and then 
describe the ways of inclusion of more details. We think that calculations in this field can be made inclusively without 
dwelling on a particular scale. Our calculations are time effective and can provide a desired accuracy. They show 
that the seemingly insurmountable complexity of geometry of the DNA in different states may be tackled because the 
geometrical differences, shown in Fig. 4, are insignificant. We would like to encourage the experimentalists to provide 
data more relevant to the actual conditions of irradiation, especially on the smallest scales involving the DNA damage. 
This information is vital for further tuning of our approach by selecting and elaboration on the most important aspects 
of the scenario. 
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